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Triangular lattice exciton model 


Daniel Gunlyck^ and Frank Tsen^ 

Naval Research Laboratory, Washington, D.C. 20375, USA 

We present a minimalistic equilateral triangular lattice model, from which we derive electron and exciton band 
structures for semiconducting transition-metal dichalcogenides. With explicit consideration of the exchange in¬ 
teraction, this model is appropriate across the spectrum from Wannier to Frenkel excitons. The single-particle 
contributions are obtained from a nearest-neighbor tight-binding model parameterized using the effective mass 
and spin-orbit coupling. The solutions to the characteristic equation, computed in direct space, are in qualita¬ 
tive agreement with first-principles calculations and highlight the inadequacy of the two-dimensional hydrogen 
model to describe the lowest-energy exciton bands. The model confirms the lack of subshell degeneracy and 
shows that the A-B exciton split depends on the electrostatic environment as well as the spin-orbit interaction. 

PACS numbers: 73.22.Dj, 73.22.Lp, 78.67.-n, 78.67.Bf 


I. INTRODUCTION 

It is well known that the electronic properties of a crys¬ 
tal depend critically on its crystal structure. Yet, the lattice 
is usually absent from descriptions of low-energy electronic 
excitations in pristine semiconducting or insulating crystals. 
The reason the lattice can often be ignored is that these exci¬ 
tations generate bound electron and hole pairs, known as ex¬ 
citons, that are typically either confined to single sites or span 
regions much larger than the relevant lattice constant. The en¬ 
ergies of the former site-confined excitons can be estimated 
by the energetics of the excited site and an “exchange energy” 
from neighboring sitesj' These excitons have been named af¬ 
ter Frenkel and can be found in e.g. molecular solids. The sec¬ 
ond type of excitons can be described by the hydrogen model 
adapted for excitons^ and are usually referred to as Wannier- 
Mott excitons and can be found in semiconductors with large 
permittivities and dispersive charge carriers. 

Lattice effects could be important, however, in the descrip¬ 
tion of excitons in semiconducting transition metal dichalco¬ 
genides such as molybdenum disulfide (MoS 2 ). These lubri¬ 
cants, investigated in the 1960s^^^ have received renewed in¬ 
terest after the isolation of individual layers^ and the demon¬ 
stration that the monolayers, unlike their bulk counterpart, are 
direct gapiS semiconductors. An estimate for the separation 
of the electron and hole in a two-dimensional (2D) semicon¬ 
ductor is the ground state radial expectation value of the 2D 
hydrogen model^iS for excitons (r)is ^ O.OQer- nm, where 
is the relative permittivity. Relative permittivities are expected 
to be quite small in 2D materialsri and for Sr ^ 3.5, the 
expectation value {r)is becomes smaller than the lattice con¬ 
stant a in these materials, which necessitates atomistic treat¬ 
ment. The exciton binding energy provides another rough 
measure of the degree to which the excitons are confined. 
Unfortunately, there is no conclusive band edge feature in 
the optical spectroscopy data, which has resulted in different 
interpretations of observed spectroscopic featuresi^ii^— Fur¬ 
thermore, band gap predictions based on density functional 
theory depends significantly on calculation details i*^'^^'^°~ — 
Photoconductivity^ and scanning tunneling spectroscopy^^ 
data suggests a binding energy of at least 0.5 eV, providing 
further indication that the exciton radius is of the order of the 



FIG. 1. (Color online) (a) Top view of the hexagonal-tiled structure 
of trigonal prismatic monolayer transition metal dichalcogenides. (b) 
The transition-metal sites form a triangular-tiled structure that ex¬ 
poses the underlying equilateral triangular lattice, (c) First Brillouin 
zone of the equilateral triangular lattice with symmetry points F, K, 
K', and M. (d) The conduction and valence band d-orbitals at K and 
K' have magnitudes that are azimuthally symmetric. 

lattice constant. 

We present a minimalistic exciton model, which we refer to 
as the triangular lattice exciton (3-ALE) model. This model 
accounts for both the exchange and spin-orbit interactions and 
the equilateral triangular lattice shown in Fig.[T] A dielectric 
constant or function is used to account for dielectric polariza¬ 
tion. Unlike typical exciton models solving the Bethe-Salpeter 
equation^S^ in reciprocal space, we derive a characteristic 
equation on sparse form in direct space, thus allowing great 
computational efficiency. This efficiency enables large calcu¬ 
lations and hence avoids the convergence challenge for tightly 
bound electron/hole pairs<^ The long-range Coulomb inter¬ 
action in direct space causes no problem for tightly bound 
electron/hole pairs in transition-metal dichalcogenides. Our 
model illustrates the breakdown of the 2D hydrogen model for 
low permittivity through the atomic scale electron-hole sep¬ 
aration and the binding energies and oscillator strengths of 
the lowest exciton states. In addition, the model shows that: 
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(1) optically allowed exciton states depend sensitively on the 
shared-site potential, as well as the relative permittivity, (2) 
the degeneracy within lower shells is broken, and (3) the en¬ 
ergy separation between the A and B excitons depends on the 
electrostatic environment, as well as the spin-orbit coupling.^ 

The 3-ALE model can, as is usually the case for exciton 
models, be separated into electron and hole contributions and 
a contribution from their mutual Coulomb interaction. Sec¬ 
tions |II] and Un] present the single-particle triangular lattice 
model and the complete 3-ALE model, respectively. Sec¬ 
tion |IV]provides conclusions about the 3-ALE model together 
with suggestions for possible generalizations. 

II. TRIANGULAR LATTICE MODEL 

Monolayers of trigonal-prismatic transition-metal dichalco- 
genides have a hexagonal-tiled structure as shown in Eig.fTfa) 
with trigonal point group symmetry D^h- The bands around 
the Eermi level are transition metal d-bands^ with only 
small contributions from the chalcogenides entering through 
hybridization^^— which we will herein neglect. This al¬ 
lows us to construct a basis comprised solely by transition- 
metal sites. Without the chalcogenide sites, we have exposed 
the equilateral triangular lattice in these crystals shown in 
Eig.HJb). This triangular-tiled structure has a point group 
pseudosymmetry D^h- Pseudosymmetry can be useful, but 
it is imperative that one exercise caution as pseudosymmetry 
introduces artificial symmetry, in the present case inversion 
symmetry. 

The triangular lattice is generated by a set of translation 
vectors R, which are linear combinations of the lattice vectors 
Si = ax and 02 = + ■\/3y)/2 with a being the lattice 

constant. We can introduce an electron field on the triangular 
lattice through the hermitian conjugate of the electron field 
operator 

^{f) = J2c„^H(l>nair- R), (1) 

naR 

where is the annihilation operator for an electron in 

the Wannier orbital^ (pna {S — R) centered on site R with 
n G {c,v} and a G {±1/2} being band and spin indices, 
respectively. Within the tight-binding approximation, these 
Wannier orbitals can be viewed as isolated atomic orbitals. 
The electronic states outside the band gap near the band edges 


TABLE I. Transition-metal dichalcogenide parameters provided in 
units of electron volts, electron mass, and nanometers. The hop¬ 
ping parameter for each material obtained from the effective mass 
m* and lattice constant a using t = 2h^/3m*a^. A is the spin-orbit 
coupling. 



m* 

a 

t 

A 

M 0 S 2 

0.64 

0.312 

0.82 

0.152 

MoSea 

0.70 

0.325 

0.69 

0.196 

WS 2 

0.49 

0.313 

1.06 

0.425 

WSe 2 

0.54 

0.325 

0.89 

0.463 



r M K r K' M r 

Wave vector 


FIG. 2. (Color online) (a) The electron band structure of the min- 
imalistic lattice model captures the band gap, band curvature, and 
valence-band split at the K and K' points, (b) Additional param¬ 
eters are needed for other band features such as a spin-split con¬ 
duction band and a valence band with a local maximum at F near 
its band edge. Bright and dark curves indicate spin. Parameters: 
Eg = 2.6 eV, t = 0.89 eV, and A = 0.463 eV. 

are concentrated around the corners of the hexagonal Brillouin 
zone,*^- at the symmetry points K and K', shown in Eig.lHc). 
These symmetry points are located at Kt = Airri/Sa with 
T G {±1} being a valley index. At these symmetry points, 
the d-orbitals in the conduction and valence bands take the 
form dz 2 and ± idxy, respectively, which are shown 

in Eig.[T}d). These orbitals are generally hybridized^i^^i^ 
throughout the Brillouin zone. This hybridization can be ig¬ 
nored, however, so long as we restrict the usage of our model 
to an energy regime near the band edges. With this restriction 
in mind, we shall herein treat the the conduction and valence 
band as two non-interacting bands. We also restrict the inter¬ 
action between different localized Wannier orbitals to nearest 
neighbors, which relative to an arbitrary site are located at the 
six sites i 5 G {±ai, ±02, ±(ai — 02)}- The Hamiltonian could 
then be expressed as H = '^^la with 

RS R 

where and £„ are nearest-neighbor and onsite hopping 
parameters, respectively. In the absence of the spin-orbit in¬ 
teraction, is spin-independent and isotropic. The latter 
isotropy is evident from the lack of azimuthal dependence of 
the magnitude of the d-orbitals in Eig.lHd). To account for 
the spin-orbit interaction, which is significant in the transition 
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FIG. 3. (Color online) Real-space probability density of the electron for a centered hole in zone-center exciton states. The triangular 

lattice grid shows that the electron and hole in the a = 1 (Is-A) and a = 2 (Is-B) exciton states are tightly bound and are likely to be found 
on the same or neighboring sites. Parameters: t — 0.89 eV, A = 0.463 eV, = 3.2 and Awq = 1.6 eV. 


metals, we add an imaginary term so that 


^naS = in+ 4icrf„ sin K+ ■ S, (3) 

where and are real parameters. We have determined 
the introduced parameters above, up to an arbitrary reference 
energy, from the band gap Eg, the effective mass m*, and the 
spin-orbit coupling A: 

Sc — “F Eg tc — t tc — 0; 

= -3f - A/2 ty = -t iy= A/18, ^ ’ 

where t = 2 Tt?/“ im*, and the band indices c and v refer to 
the the conduction and valence band, respectively. 

The single-particle band structure 

+ X! ^naS (5) 

S 


for an electron and hole, respectively, we can express a general 
excitation state as 

lexc) = —^ f - |0)) (7) 

' ' ^^crKXR c!R'+R SR'' '' ' ^ 

RR’ 


where N is the number of lattice sites, direct 

space eigenfunction for the electron/hole pair, and |0) is the 
vacuum state. The focus below is to obtain i^g^xR solving 
for the relative motion of the electron and hole. 

The Hamiltonian H describing the electron-hole pair con¬ 
tains a mutual Coulomb interaction as well as the single¬ 
particle contributions described above. We neglect Coulomb 
contributions from electrons in the valence bands and approx¬ 
imate the screened Coulomb integrals as 


Vs=l ^-^reo\R\ 

1 Avo iR = 0), 


( 8 ) 


is shown in Fig. |2] for tungsten diselenide (WSe 2 ). Parame¬ 
ters for additional semiconducting transition metal dichalco- 
genides are provided in Table |I] 


III. TRIANGULAR LATTICE EXCITON MODEL 


When dealing with elementary excitations, it is convenient 
to perform a canonical transformation that takes into account 
electrons in the conduction band and holes in the valence 
band. Consider an elementary excitation bringing an electron 
with spin a from a state in the valence band with wave vector 
ky = k — K/2 across the band gap into a state in the conduc¬ 
tion band with wave vector ky = k + K/2. This excitation 
creates an electron with fee = kc and a hole with kh = —ky 
so that the combined wave vector becomes K. This combined 
wave vector commutes with the Hamiltonian and is therefore 
a good quantum number. Other quantum numbers are a and 
A, the latter describing the relative motion of the electron and 
hole. Dehning the phase-modihed annihilation operators 


d-s = crc' -e 

vaR 


-iK-R/2 

5 

-iK-R!2 

5 


(6a) 

(6b) 


where e is the elementary charge, Sq is the vacuum permittiv¬ 
ity, and Auq = vq — wq is the difference between the direct 
integral vq, and the exchange integral wq, the latter assumed 
to be negligible when the electron and hole occupy different 
sites. 

In analogy with the hydrogen atom, we choose the origin to 
be position of the positive charge carrier, in this case, the hole. 
We also let this origin track the motion of the hole, thereby 
transfering the kinetic energy of the hole over to the electron. 
This allows us to express the hole-centered electron/hole pair 
Hamiltonian as 

^crii — X/ '^(TKS^\g+^aR + '^(^0 - ^R)c\gC^g, (9) 
m R 


where the coefficients 

^aKS — ^caS 

Eq = Sy £y Eg , 


— t ?e 

verb 


iK-5l2 


(lOa) 

(10b) 


can be obtained from Eqs. Q and (|4]i with parameters from 
TableU 

The solutions of the characteristic equation for the relative 
motion 


H^g\aKX)=E^g^\aKX) (11) 
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FIG. 4. (Color online) Zone-center excitons. Exciton energy as a 
function of the relative permittivity (a) and onsite Coulomb interac¬ 
tion (b) with the intersection Sr ~ 3.2 and Ahq = 1.6 indicated 
by vertical lines. Solid and dashed curves obtained from the 3-ALE 
model and the the 2D hydrogen model, respectively. While there is 
good agreement for the a = 3,4 (2p) curve, the latter model breaks 
down for the tightly bound electron/hole pairs in the a = 1 (Is- 
A) and a = 2 (Is-B) states. The energy separation between these 
states depends on the relative permittivity (c) and onsite Coulomb 
interaction (d), as well as the spin-orbit coupling A. The oscilla¬ 
tor strengths relative to that of a = 1 (Is-A) are shown in (e) and 
(f). Dark and bright curves represent exciton states from the A series 
and B series, respectively. Additional parameters: t = 0.89 eV and 
A = 0.463 eV. 


provide the electron/hole pair eigenenergies and eigen¬ 

states of the form 

( 12 ) 

R 

These states can be divided into two sets: (1) bound 
electron/hole-pair states known as excitons and (2) unbound 
electron/hole-pair states. The former states appear at discrete 
energies and are described by A = a, where we let a take pos¬ 
itive integers counted from the strongest bound state. The lat- 


TABLE II. Hydrogen shell labels for the lowest five quantum num¬ 
bers describing excitations from the valence band (Series A) and 
split-off band (Series B) resulting from the spin-orbit coupling. 


Quantum number 

a = 1 

a = 2 

a = 3 

Q = 4 

a = 5 

Series A label 

Is 



2py 

2s 

Series B label 


Is 





0.25 


0 

> 

A. 

^ -0.25 

0) 
c 
LU 

-0.5 

-0.75 

r M K r K' M r 

Wave vector 



FIG. 5. (Color online) Electron/hole pair band structure. Discrete 
bands and filled regions contain bound exciton states and free elec¬ 
tron/hole pairs, respectively. Bright and dark bands indicate spin of 
the excited electron. The bands, enumerated at F starting from the 
strongest bound state, are: a = 1 (Is-A), a = 2 (Is-B), a = 3, 4 
(2p), and a = 5 (2s). Parameters: t = 0.89 eV, A = 0.463 eV, 
Et = 3.2, and Adq = 1.6 eV. 


ter unbound states are continuous in energy and are described 
by A = fc. 

The exciton states with the wave vector K at the zone center 
r are of particular importance for optoelectronic experiments 
due to the relative high speed of light causing almost vertical 
electronic excitations. The eigenfunctions of the hve strongest 
bound r excitons are shown in Fig. [3 where each vertex rep¬ 
resents a transition metal lattice site. Strictly, the lattice sym¬ 
metry prohibits the exciton quantum numbers a from being 
separated into radial and azimuthal quantum numbers. For the 
lowest-energy bands, however, designation using shell labels 
from the hydrogen model still makes sense and we provide the 
conversion in Table HJ 

The exciton energies depend on the screened Coulomb in¬ 
teraction, as depicted in Fig.|4] First, we note that the a = 3,4 
(2p) curves are lower in energy than the a = 5 (2s) curve, in 
agreement with hrst-principles predictions jF® Second, we note 
that the a = 2 (Is-B) exciton only appears below the band 
edge for a sufficiently small dielectric constant. In this regime, 
the hydrogen model also breaks down for the A exciton series 
due to the tightly bound electron/hole pairs. Third, the a = 
3,4 (2p) states have a node at the origin and are unaffected by 
Aug. Fourth, the energy separation between the a = 1 and 
a = 2 states (“the well-known A-B exciton split”) is not only 
a function of the spin-orbit coupling but also the Coulomb 
interaction. Fifth, the Coulomb interaction also affects the 
relative oscillator strengths A^jAx = ICoaol^/ICoiol^ 
could be modulated through the dielectric environment. 

Exciton states away from the zone center could also po¬ 
tentially be exploited in applications. The energies of these 
states, E^g^, are shown as discrete bands in Fig.|3 Reflection 
symmetry or time-reversal symmetry causes the zone-center 
states to be doubly degenerate and requires the K excitons to 
have the same energies as the K' excitons with the opposite 
excited spini^2i2^ This degeneracy could be broken, however, 
by an external magnetic field42ili Transitions between exciton 
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states are possible through additional interactions, as long as 
they exhibit appropriate symmetryFirst-principles band 
structures indicate that there might be additional states near 
K and K' originating from the valence band F point. To cap¬ 
ture these states, the single-particle basis would need to be 
expanded. 

The electron/hole pair band structure also contains unbound 
states. Unbound states have been derived analytically for the 
2D hydrogen model.Obtaining similar solutions using the 
3-ALE model is not practical, owing to the large number of 
sites needed for such calculations. Therefore, we model in¬ 
stead the filled regions in Fig.|5]using free electron/hole pairs. 
For these free electron/hole pairs, the Coulomb interaction 
^ 0. The eigenfunctions are then Bloch waves given by 

^crKkR ~ The corresponding energies of the 

free electron/hole pairs are 

S 

where we have as expected = E^_g_^, again from 

time-reversal symmetry. 

IV. CONCLUSIONS 

The model presented herein have been derived from the 
many-body Hamiltonian with the goal to extract the dominant 
physics among the elementary excitations in 2D semiconduc¬ 
tors with a equilateral triangular lattice. The approximations 


should be reasonable for tightly bound electron/hole pairs, 
and we therefore expect the model to work well for optically 
driven transitions into the lower exciton states. The predicted 
results are in qualitative agreement with first-principles quasi¬ 
particle calculations, and for improved accuracy, we suggest 
using this model with a more detailed dielectric environment 
with a complex dielectric function that is spatially^ and/or 
energy dependent4^ To achieve a better description away 
from the band edges, which could be needed for zone-corner 
excitons, it is also recommended to use a larger basis to ac¬ 
count for d-band hybridizationi^iiii^ During the final stage 
of this paper, the authors have become aware of a recent paper 
including such hybridization effects for excitons in MoS 2 ^ 
In summary, we have presented a minimalistic and com¬ 
putationally efficient model capturing the essential physics 
of low-energy electronic excitations. This model could be 
used to understand experimental absorption spectra and is 
well suited to serve as a foundation for more advanced mod¬ 
els, which could describe exciton coupling to, e.g., defects, 
interfaces, vibrations, and externally applied fields. 
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